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We evidence the importance of electron charging under nonequilibrium conditions for carbon- 
nanotube-based molecular bridges, using a self-consistent Green's function method with an extended 
Hiickel Hamiltonian and a three-dimensional Poisson solver. Our analysis demonstrates that such 
feature is highly dependent on the chirality of the carbon nanotube as well as on the type of the 
contact metal, conditioning in a nongeneralized way the system's conduction mechanism. Based on 
its impact on transport, we argue that self-consistency is essential for the current- voltage calculations 
of semiconducting nanotubes, whereas less significant in the case of metallic ones. 



The fabrication of carbon nanotube (CNT) field-effect 
transistors [JJ has seen a constantly increasing trend in 
the last years, due to the advancements in the character- 
ization and manipulation techniques, as well as the in- 
dustry's ongoing request for prototype devices that could 
allow for a shift from the actual silicon-based to that 
of carbon nanoscale electronics. Alongside, the need 
to theoretically interpret the behavior of such compo- 
nents has attracted a large amount of research, leading 
to an optimization of the production process in the lab- 
oratory. It is therefore well-known by now that equilib- 
rium electron charging effects can play an important role 
in the operation of CNT-based molecular bridges, pro- 
voking Schottky-type potential barriers of non-negligible 
dimensions (21 13] in the interface area between the carbon 
nanotube and the contacts. Electron charging however is 
also present under non-equilibrium conditions [3] due to 
both spatially and energetically anisotropic interactions 
between the source and drain electrodes and the CNT's 
local density of states at energies close to the electro- 
chemical potentials of the contacts. This type of phe- 
nomenon, nonetheless significant for the quantum trans- 
port mechanism of systems within the nanoscale, has re- 
ceived a minor attention in the carbon nanotube context. 

In order to examine charging phenomena between 
CNTs and real metallic leads, approaches that go beyond 
single-electron tight binding models need to be adopted. 
In this sense, ab initio models have established accu- 
rate system descriptions j5j [6j [7], bearing though a non- 
negligible computational load. Alternatively, semicm- 
pirical extended Hiickel approximations have demon- 
strated various merits when employed for the extrac- 
tion of CNT transport attributes (e.g. secondary en- 
ergy gaps in zigzag metallic nanotubes [8], feasible nar- 
row diameter CNT study [9] etc.), while being capable of 
a reliable description of molecular conduction under non- 
equilibrium^ in a more affordable way. Here we present 
self-consistent quantum transport calculations for device 
structures based on finite semiconducting and metallic 
CNTs. We couple the non-equilibrium Green's function 
formalism [10] (based on an extended Hiickel Hamilto- 
nian) with a full 3D Poisson solver for a realistic repre- 



sentation of the device and contacts' chemistry and the 
system's electrostatics. The goal of this study is to evi- 
dence the role of electron charging under nonequilibrium 
conditions in the presence of different contact metals (Au, 
Pt and Al), and focus on the variations of the system's 
transmission probability when tuning the terminal po- 
tentials, which represent the external parameters of our 
model. 

Our approach is based on the single particle retarded 
Green's function matrix G = [ES — H — — 
where E is the scalar energy, H is the 'device' Hamilto- 
nian matrix in an appropriate basis set, S is the overlap 
matrix in that basis set and ^l,r is the self energy, which 
includes the effect of scattering due to the left (L) and 
right {R) contacts. A Landauer-type expression can be 
used for the current calculation in case of coherent trans- 
port: 
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I= ^J dET(E)[f(E,» L )-f(E,» R )], (1) 

where T(E) = Tr[T lGT rG'] is the transmission as a 
function of energy, T L>R = i[^L,R-^ LtR ] and f(E, p-l.r) 
represents the Fermi-Dirac distribution of electrons in 
the contact at chemical potential pl,r (T=300X in this 
work). For the description of both device and contacts 
we use an Extended Hiickel semiempirical Hamiltonian 
calculated in a non-orthogonal basis set of Slater-type 
orbital functions [9]. Charging effects can be introduced 
in the formalism with the inclusion of a self-consistent 
potential Usc(Ap) that is added in the bare device's 
Hamilonian H (H = H + U sc {Ap)). The U SC {Ap) 
term can then be determined by the approach of Zahid 
et aZ[4j: 

Usc(Ap) = U L aplace + Up 

oisson 

(Ap) + U Image (Ap), (2) 

where Ap represents the change in the charge density 
between the nonequilibrium and the equilibrium condi- 
tions (Ap = p — p eq ), and p is given by the expression 
below: 
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FIG. 1: Geometry, potential profile and change of the elec- 
tronic density Ap for a 2-unit cell (3,2) CNT with Au(ffl) 
metallic contacts, when a 2V source-drain bias is applied. 
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P=— dE[f(E, fi L )GT L G^ + f(E, p R )GT R G^} 
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(3) 

The calculation of the Poisson term can derive in the 
framework of the complete neglect of differential orbital 
theory, using only the Hartree potential for the Coulomb 
interaction [4 . Laplace and Image expressions are deter- 
mined numerically by solving the V 2 C7 = equation with 
a finite element method in real space, using the appropri- 
ate boundary conditions |4J. All three potential compo- 
nents are evaluated on the atomic sites of the CNT. Fi- 
nally, equation [3] results computationally demanding and 
therefore contour integration techniques in the complex 
energy plane [TX] for energies smaller than min(fiL, p R ), 
as well as Gaussian quadrature formula implementations 
have been introduced in the model for optimization pur- 
poses. 

We have considered molecular bridges based on defect- 
free single-wall carbon nanotubes and restricted our cal- 
culations in the ballistic regime, since the maximum 
length of the device tubes is not greater than a few 
nanometers, less than mean free path measurements [L?j 
even for high biases. As the focus of our investigation 
lies in the nonequilibrium regime, we have separated the 
noncquilibrium transport features from the electrostatic 
effects of equilibrium (e.g. Schottky barriers) by impos- 
ing the alignment between the metal's work function and 
the charge neutrality level of the CNT for both equilib- 
rium and noncquilibrium electronic density calculations. 
This condition can be practically realized shifting the 
CNT energy bands with respect to the metallic ones, 
similarly to using a gate electrode to control a CNT field- 
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FIG. 2: Potential profile of a (3,2) CNT when a IV (a) and 
3V (6) source-drain bias is applied, for Au(lfl), Al(IIl) and 
Pt(lll) metallic contacts. 



effect transistor. We have, moreover, set the distance of 
the source and drain electrodes from the CNT ends to 
lA in order to ensure a device functionality in the self- 
consistent field regime and avoid undesirable weak cou- 
pling effects induced by the finite size [5], which exceed 
the current study's objectives. 

Figure [T] shows the composite potential profile of a 
(3,2) CNT with Au(lll) contacts when a 2V bias is ap- 
plied, as well as the change of the electronic density Ap 
compared to equilibrium conditions. The potential terms 
that arise due to nonequilibrium charging are the Pois- 
son and Image ones, which consequently compromise the 
overall term, correcting it with regard to the electrostatic 
Laplace level. Oscillations of the potential calculated 
at the atomic sites are due to a higher accumulation of 
nonuniformly distributed charges in the interface regions 
that reflect the axially nonsymmetrical interface atom 
positioning. Our analysis has demonstrated that there 
are two main factors that shape the aforementioned po- 
tential components: a) the chirality of the carbon nan- 
otube and b) the type of the metallic contacts (figure 
[2| . A source-drain bias quantitative dependence has also 
been observed (see figure ^p) . It is interesting to note 
that the potential drop in the interface area between the 
nanotube and the contact is also element-dependent and 
that charging can break the profile's symmetry with re- 
spect to the center of the CNT, as can be clearly seen 
in the case of the 2>V bias. Figure [3] shows how nonequi- 
librium charging affects the transmission probability of a 
(3,2) and a (9,0) system in the presence of Au leads. The 
transmission of the respective bulk systems has also been 
plotted in order to clarify that nonetheless the small size, 
the evanescent modes of the system [IB] do not decisively 
determine the transport features with respect to the long 
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FIG. 3: Transmission as a function of energy for a (3,2) 2-unit 
cell (left column) and a (9,0) 5-unit cell CNT (right column), 
for 0.2V, IV and 3V source-drain biases. The contact metal 
is Au(lll) whereas zero energy refers to the charge neutrality 
level of each CNT. Interrupted lines represent the transmis- 
sion probability of respective CNTs in the bulk limit under 
equilibrium. 



nanotubc limit. In both cases and as the bias rises up 
to 3V, we can observe a progressive spatial redistribu- 
tion of the transmission peaks throughout the energy 
spectrum, which is a direct consequence of the alteration 
of the system's energy eigenvalues. Although such fea- 
ture is common for the semiconducting and the metallic 
CNT, their impact on the conduction mechanism is evi- 
dently distinct. Whereas in the metallic tube's case the 
alteration of the transmission curve does not change the 
principal transport characteristics, in the semiconduct- 
ing case we can observe a clear unidirectional motion of 
the valence band towards smaller energies, which results 
in a widening of the conduction gap (from about 0.88ey 
with 200mV bias to lAAeV with 3V bias). Analogous 
responses are also obtained with the other contact met- 
als, although the positioning of the conduction band is 
moreover contact dependent. 

Charging can be perceived as an interaction between 
the source and drain electrochemical potentials and the 
device's local density of states. Since the latter is differ- 
ent in the valence and conduction band zones (see the 
transmission functions of figure [3j, the contacts' elec- 
trochemical potentials preferentially correlate with the 
states of the one band rather than the other, provok- 
ing a charge accumulation on the CNT body that can 
be only captured by a self-consistent approach. This 
has a practical consequence on the characteristics of the 
calculated current-voltage curves of the studied CNTs. 
As we can see in figure [4] in the semiconducting tube's 
case, a substantial reduction of the overall source-drain 
current can be noticed with respect to the calculated 
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FIG. 4: Current-voltage curves for a (3,2) 2-unit cell in the 
presence of Au(lll) (a), Pt(lll) (6) and Al(lll) (c) contacts, 
and for a (9,0) 5-unit cell CNT in the presence of Au(lll) 
contacts (d). Dotted lines represent the respective I-V curves 
calculated without self-consistency. 

values without the inclusion of self-consistency (based 
on the simplified potential model of ref.[9l [H]). Even 
if the effect can be observed for all contact metals, its 
impact is variable, being higher in the Pt case rather 
than in the Au and Al ones. On the contrary, in the 
metallic tube's case, charging has a minimal influence on 
the overall conduction mechanism, making possible valid 
transport calculations without the need to include self- 
consistency. Finally, early calculations (not shown) for 
semiconducting and metallic CNTs of different helicities 
demonstrate a qualitative correspondence to the afore- 
mentioned scheme, although results are quantitatively 
distinguishable from case to case. 

To summarize, we have seen how nonequilibrium 
charging can influence the conduction mechanism of 
CNT-based molecular bridges, demonstrating that the 
latter depends on both the chirality of the CNT as well 
as on the type of the metallic contact. We have moreover 
showed that the corrections induced in the calculation of 
the I-V curve are essential in the case of the semicon- 
ducting tube as well as significantly contact-dependent, 
whereas less significant in the metallic tube's case. The 
importance of such findings can reflect on the common 
perception of transparency for metals when used as con- 
tacts to carbon nanotubes. Our analysis has demon- 
strated that this cannot be defined unilaterally by the 
height of the Schottky barrier formation, if other contact- 
induced features that have an impact on transport are 
not evaluated. Moreover, apart from carrier-movement 
phenomena, quantum-chemical effects related to the cou- 
pling between the CNT and the metallic electrode(e.g. 
blocking of conduction channels [9 ) can influence strongly 
the current-carrying capacity of a carbon nanotubc. The 
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extent of all aforementioned characteristics is intrinsic 
to the CNT type. This means that no generalized con- 
clusions can be drawn regarding the quality of a metal- 
lic contact when employed in such context. Rather, a 
more feasible categorization process could be based on 
self-consistent results obtained for CNTs of similar elec- 
trical and geometrical properties (conducting character, 
charge neutrality level, diameter). 
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